Controlling growth of nanorod via screening effects 
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We demonstrate in this report that nanorods can spontaneously grow on the top of a transient 
mound when the screening effect is considered in island growth. The number of the topmost growing 
layers is defined as the active length, which decreases as the screening effect is enhanced. It follows 
that the thickness of the transient mound underneath the nanorod increases as the square of the 
active length, and a smaller active length favors the formation of nanorod. 

PACS numbers: 61.46.-w, 68.65.Ac, 68.55.-a 
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Due to the size confinement of nanorods, their unique 
electric and optical properties make them prospective 
materials for device applications in nanoelectronics and 
nanophotonics [ll-Q. However, controlling the growth 
of nanorod is usually tricky, and the understanding of 
the associated growth mechanism remains a challenging 
task. In metallic growth system, for example, in the ab- 
sence of catalysts, instead of generating nanorods, very 
often a structure calld "wedding-cake" or "mound" is 
generated (3, 13- To solve this problem, oblique angle 
deposition is normally used, in which the upmost sur- 
face features shadow nearby lower surface regions where 
the incident atoms cannot reach In addition, in 

some systems the long-range attractive forces between 
the incoming atoms and grown structure inherently steer 
the incident atoms to the upper growth sites, and hence 
stimulate the growth of nanorods in these regions [13] ■ 
The essential feature in these two scenarios is that the 
topmost structures are superior in capturing incoming 
atoms. Here we term such a phenomenon as the screen- 
ing effect. In fact, it is a well established effect in Lapla- 
cian growth: Among the protruding growing parts, the 
front-most tips get more nutrients and hence grow faster. 
Although screening effect has been observed to suppress 
the growth of those being trapped behind and thus form 
more ordered filament structures 11, l^], the quantita- 



fiuence of surface kinetics on the growth mode and radius 
selection of nanorod growth 15 j . 



tive understanding of the screening effect in the inter- 
facial growth, especially the role of screening effect in 
keeping the geometrical shape of the nanorods, has not 
been well addressed so far to the best of our knowledge. 

In general, interfacial growth involves landing and 
accommodation of atoms on the surface, diffusion of 
adatoms on the surfaces, nucleation, and formation of a 
new crystalline layer by lateral expansion of the nucleus. 
The screening effect influences the landing positions of 
the deposit atoms, and surface diffusion promotes their 
redistribution. So both screening effect and surface ki- 
netics govern the distribution of the adatoms on crystal 
surface, and affect the nucleation and growth processes. 
The influences of surface kinetics on the growth mor- 
phology and size of nanostructures have been extensively 
studied 

[iMIIl- Recently we also demonstrated the in- 



In this report we focus on the influence of screening 
effect on the formation of nanorods with uniform cross- 
section size. Our calculations indicate that the nanorods 
can spontaneously grow on top of a transient mound 
when the screening effect has been considered in island 
growth. The number of the topmost growing layers is de- 
fined as the active length, which decreases as the screen- 
ing effect is enhanced. Consequently, the thickness of the 
transient mound underneath the nanorod increases as the 
square of the active length, and a smaller active length 
favors the formation of nanorods. 

Consider a growing stepped nanostructur with n 
atomic layers. The dimensionless area of the substrate 
occupied by the nanostructure is denoted as Aq, and that 
of the i-th layer relative to the substrate is as Ai . In is- 
land growth, the energy barrier at the step is usually 
larger than that within the terrace. The growth of a 
partially buried layer is maintained by the directly de- 
posit atoms in the underneath layer, since the interlayer 
transport of the adatoms on the buried layers can be neg- 
ligible. Assuming adatoms are deposited with rate in the 
normal direction of F per surface lattice site. If there 
is no screening effect, the deposition is uniform on the 
surface, and all the layers can keep growing by captur- 
ing the deposit atoms. The opposite limiting case is that 
only the topmost layer is allowed to grow as a result of 
extremely strong screening effects. The more realistic sit- 
uation is that the moderate strong screening effect leads 
to a finite number of the active layers, i.e., the topmost 
layers which are active in capturing the deposit atoms. 
Correspondingly, only the topmost Ng layers keep grow 
as the number of atomic layers n adds up with deposi- 
tion. Therefore the quantity Ng can be regarded as a 
measure of the screening strength: The smaller Ng is, 
the stronger the screening is. We refer hereafter Ng as 
the active length. 

By defining the dimensionless time as T = Ft at depo- 
sition time t, the coverage 9i, which is defined as Ai/A^, 




FIG. 1: Morphologies without screening effect {Ng = oo, up- 
per panels) and Ng — 20 (lower panels). (a,c) Surface profiles 
of morphologies at different T, with characteristic coverage 
dc = 0.1. (b) Convergence of the normalized profiles when 
without screening effect, depending on different 9c. (d) Nor- 
malized profiles with Ng = 20 and 9c = 0.1 at different T. 



increases with time according to the rate equation [l^ , 

de. 



dT 
d0^ 
dT 



= it<n~Ng), (1) 
= e,_i-0, (n-Ng<i<n). (2) 



Since there are no sinks for atoms atop, the topmost ter- 
race absorbs all the atoms landing on layers n and n — 1 , 
and grows according to 



dK 
dT 



— On-l- 



(3) 



The influences of interlayer hoppting rate u' and the 
deposition rate F are included into a characteristic cov- 
erage which is defined as 9c = Ac/Aq, where Ac = 
{ly/TTv' /F)^^^ is the dimensionless characteristic area ac- 



cording to Ref. [15[. Starting from the substrate with 



area Aq larger than Ac, the coverage of the i-th layer 6i 
at different time can be obtained by integrating Eq. ^ 
with a time step AT. In this work, we use AT = 0.01, 
and the influence of different AT is discussed later. 

In the general situation of growth, the time interval be- 
tween subsequent deposition events is much larger than 
the survival time of the adatom on A„. For the simplest 
situation where the smallest stable nucleus is a dimer, if 
a circular cross section is assumed, the nucleation rate on 

top of An per unit time can be deduced 51 = 2^^' Ml ■ 
The corresponding nucleation rate per unit dimensionless 
time is fir = ^/F = 1{A„/Ac)^/^. The random nucle- 
ation process is then modeled by introducing a uniform 



distributed random number between and unit: When- 
ever fJ^AT is larger than the random number, a new 
layer is added atop and n increases by one. 

The variation of 9i is recorded under specific growth 
conditions. The surface profiles without screening effect 
at different T are shown in Fig. 1(a). The lateral position 
on the X axis is proportional to It shows that the 

number of the exposed layers increases with time. After 
re-scaling, the height profile of the mound converges to a 
time-independent asymptotic shape which only depends 
on 9c, as shown in Fig. [Ijb). This is a typical character- 
istic of the morphology known as wedding-cake (jl. [l6l|. 

The general evolvement of the surface profile when 
screening effect is introduced is shown in Fig. [Ijc) , where 
Ng — 20 and = 0.1. It indicates that the coverage de- 
creases rapidly as a function of i, until it approaches 9c 
within nt transient layers. The corresponding re-scaled 
profiles at different T are given in Fig. [TJd). It shows 
that the curve of T < Ng is superposed on the the one 
without screening effect, which suggests that the island 
grows just as in the situation of infinite Ng, keeping the 
well-known morphology of wedding-cake shape. Once n 
is larger than iV^, it is obvious that the re-scaled profiles 
in Fig. [ifd) gradually depart from the wedding-cake one, 
since the bottom n — Ng layers cease growing. 

To identify the physical meaning of Ac, we repeat 
the numerical simulations with different 9c while keep- 
ing Ng a constant. For the stable morphologies ob- 
tained after sufficient growth, we find that \ii{\n{0i/9c)) 
decreases approximately linearly as a function of i, as 
shown in Fig. EJa), until it decreases to a constant C 
when i — nt, and then it remains. It suggests that a 
nanorod with area of A^. = exp(exp(C)) is grown on 
top of the mound. The constant C decreases with de- 
creasing AT used in the simulations. When AT = 0.01, 
C ~ — 5 which corresponds to a stable nanorod with area 
Ar ~ 1.00676Ac. It justifies the importance of Ac as a 
characteristic area [15|. Realistically, Ar should always 
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FIG. 3: Logarithmic plot of the slope of the linear region of 
the curves shown in Fig. [2] with 9c = 0.01 (dots) and that of 
K = ln(3.5)/Afg (dashed line). 



be integral in unit of a§, therefore Ai/Ac decreases and 
approaches to Ar/Ac ~ (1 + as the limit. For the 

Ac in nanoscale, C is in the range of —4.5 ^ —6.9. 

In order to study the quantitative influence of the 
screening effect, the variation of coverage 9i for differ- 
ent Ng is simulated by setting 0c = 0.05. The re-scaled 
coverage \n{\ii{9.i/0c)) are shown in Fig. [2Jb), as a func- 
tion of i. The main features are similar as that discussed 
above for Ng = 20 in Fig. [2ja). The screening effect 
just changes the magnitude of the linear slope: When 
Ng is increased as a result of weaker screening effect, the 
magnitude of the slope becomes smaller correspondingly. 

The dependence of slopes on Ng can be obtained by 
fitting the linear region of the curves in Fig. [Ifb) . The 
results are plotted logarithmically in Fig. [3l where K is 
positive denoting the magnitude of the slope. It reveals 
that K is inversely proportional to iV^ , 



K = ln(3.5)/7V, 



The area of Aj, is then obtained as, 

Ac 9c ' 



{rii < i < nt) 



(4) 



(5) 



The pow index becomes (|)' when Ng = 1, which is 
consistent with the analytic results in Ref. |15| . 

It is thus clear that starting from the substrate with 
cross sectional area Aq larger than Ac, a nanorod with 
area of Ar is finally grown on top of a transient mound 
of nt layers. According to Eq. (I5|), the number of the 
transient layers nt is proportional to N^ , with the coeffi- 
cient a ~ [-\n{\TL{Ar/Ac)) - ln(ln(AoMc))]/ ln3.5. The 



value of a is in the range of 2 '-^ 7 for 9c from 0.001 to 
0.9 and Ac in nanoscale. The growth of a nanostruc- 
ture therefore shows different feathers at different du- 
ration: Wedding-cake morphologies {n < Ng), tapered 
morphologies {Ng < n < aNg) and nanorods with a ta- 
pered base [n > aNg) can be observed successively. Ac 
determines the area of the approached nanorod, whereas 
the screening strength determines the transient layers, 
nt — aNg. At the limiting case when no screening exists, 
the well-known wedding cake morphology is obtained. 

In conclusion, we have conducted simulations to study 
the screening effect on the nanorod growth. We found 
that a specific growth system can be well characterized as 
the characteristic area of Ac and the active length Ng . Ac 
determines the area of the approached nanorod, whereas 
the screening strength influences the transition process 
from taper-like morphology to the nanorod, nt = aNg. 
Larger screening effect leads to smaller active length Ng 
and favors nanorod growth. Therefore it is possible to 
control the nanorod growth process by tuning Ng via 
varying the screening strength, for example by varying 
the deposition angle in oblique angle deposition. 
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